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ABSTRACT 

We analyze the standardizability of Type la supernovae (SNe la) in the 
near-infrared (NIR) by investigating the correlation between observed peak NIR 
(YJH) absolute magnitude and post-maximum S-band decline-rate [Amis (£>)]. 
A sample of 27 low-redshift SNe la with well-observed NIR light-curves observed 
by the Carnegie Supernova Project (CSP) between 2004 and 2007 is used. All 
27 objects have pre-maximum coverage in optical bands, with a subset of 13 
having pre-maximum NIR observations as well; coverage of the other 14 begins 
shortly after NIR maximum brightness. We describe the methods used to derive 
light-curve parameters (absolute peak magnitudes and decline-rates) from both 
spline and template fitting procedures, and confirm prior findings that fitting 
templates to SNe la light-curves in the NIR is problematic due to the diversity 
of post-maximum behavior of objects that are characterized by similar Am 15 (B) 
values, especially at high decline-rates. Nevertheless, we show that NIR light- 
curves can be reasonably fit with a template, especially if the observations begin 
within 5 days after NIR maximum. SNe la appear to be better "standardizable 
candles" in the NIR bands than in the optical bands. For the subset of 13 objects 
in our data set that excludes the highly reddened and fast-declining SNe la, and 
includes only those objects for which NIR observations began prior to 5 days 
after maximum light, we find modest (1-7 a) evidence for a peak luminosity vs. 
decline-rate relation in Y, and stronger evidence (2.8a) in J and H. Using R v 
values differing from the canonical value (Ry = 3.1) is shown to have little effect 
on the results. A Hubble diagram is presented for the NIR bands and the B- 
band. The resulting scatter for the combined NIR bands is 0.13 mag, while the 
S-band produces a scatter of 0.22 mag. Finally, we find evidence for a bimodal 
distribution in the NIR absolute magnitudes of fast-declining SNe la [Ami 5 (B) > 
1.7]. These data suggest that applying a correction to SNe la peak luminosities 
for decline-rate is likely to be beneficial in the J- and if-bands to make SNe la 
more precise distance indicators, but of only marginal importance in the Y~-band. 

Subject headings: Supernovae 
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Introduction 



It is widely accepted that SNe la are excellent standardizable candles at optical wave- 
lengths. After applying an empirical correction established between light-curve shape and 
peak m agnitude , SNe la become one of the most precise extragalactic distance indicators 
known. iPhillipsI ( 119931 ) was the first to discover a tight correlation (hereafter the "Phillips 
method") between optical absolute magnitudes at maximum light and the decline-rate pa- 
rameter Amis (5), defined as the drop in 5-band brightness from peak to 15 days later. 
Slow-declining SNe la [i.e., lower Ami^B)] are intrinsically brighter than their fast-declining 
count erparts. Since Phil l ips' initial work, the Phillips method has continued to be uti- 
lized faamuv et aD Il996t IPhillips et aD Il999l : IPrieto et all l2006[ burns et al.l boilh. alon g 



with other similar cal i bratio n methods, such as the "stretch" method (IPerlmutter et al. 



1997 



Goldhaber et al. 



Riess et al. 



(CMAGIC ; 



1996 



1998 



20011). the mult icolor light-curve shape method (MLCS/MLCS2k2; 



Jha et al.l 120071 ). the color-magnitude intercept calibration method 



Wang et al.ll20031) . the s pectral adaptive lig ht-curve template method (SALT/SALT2; 



Guy et aU2005l . 120071 ). and SiFTO (jConley et al.ll2008[ ). Applying such correct ions decreases 



the d i spersion in measu r ed distance estimates to the 0.20 mag level o r less (IHamuy et al. 



19961 : IPrieto et al.l 1200ft IWood-Vasev et al.l 120081 : iFolatelli et al.ll2010h . which significantly 
improves the accuracy of SNe la as standardizable candles at optical wavelengths. 

SNe la have played a critical role in cosmological studies over the past two decades. 
While the majority of SNe la studies have been performed using optical bands, there has 
long been the hope that they might actually have greater accuracy in NIR bands, due to 
the reduced effe cts of dust extinction, which is one of the main sources of error in distance 
determinations ( jFreedman et al.l 120091 ). This has motivated recent work on establishing the 
standardizability of SNe la in the NIR, which has been shown to be advantageous for several 
reasons. Extin ction corrections are smaller, on the order of a magnitude less, in NIR than in 
optical bands ( jCardelli et al.lll989[ ). and there appears to be a shallower dependence of abso- 
lute NIR magnitude on decline-rate compar e d with optical band s (Krisciun as et al. 2004a,c; 



Wood-Vasev et al. 



2008 



Sullivan et al. 



Mandel et al. 



6b7o|), 



and 



20ld: liblatglh et alJlioiol ). Recently. iGallaeher et al. 
Kelly et al. ni- 

tudes of SNe la depend on host-galaxy properties (i.e., star formation rate, host-galaxy mass, 
host-galaxy metallicity), making the correction factors evolve with redshift. Therefore, min- 
imizing these correction factors by observing in the NIR is important, not only to decrease 
the random errors incurred for each SN la's correction, but also to reduce our sensitivity to 
any systematic evolution of the correction factors with redshift. 



The first NIR observation of a SN la was reported by iKirshner et al.l (119731 ). and the 
first extensive NIR data set was obtained by Elias et al. (1981, 1985). Since that time, there 
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has been a large increase in the number o f SNe la ligh t -curv e s observed at NIR wavelengths, 



espec ially within the last decade (e.g., iFrogel et al 



2008 



1987 



Jha et al.l 



2003 ; Krisciunas et al. 2004a,b,c, 2006; IStritzinger k, Sollerman 



Contreras et al.l l2010t IStritzinger et al.l 1201 ll and references therein). Two groups, in 



2007 



1999; 



Phil lips et al. 



Wood-Vasey et al. 



particular, have been working to obtain homogenous opt ical and NIR data sets of SNe of 
all types in multiple filters: the C fA Supernova Grou p (jWood-Vasey et al.l 120081 ) and the 
Carnegie Supernova Project (CSP; lHamuy et al.ll2006l ). 



Initial studies suggested that SNe la may present a weaker peak luminosity vs. decline- 
rate dependence in the NIR compared with the optical, which display slopes of 0.63, 0.61, 
0.57, and 0.52 ma g unit -1 decline-rate [as measured by A misf-B)] in BVRI, respectively 



(Prietoeta 



2006). This led some early investigators (e.g.. lElias et al.lll985l ; |Phillipslll993 



Meikld 120001 ) to suggest that these objects may be nearly perfect standard candles in the 
NIR. Krisciunas et al. (2004a, c) confirmed these suggestions, finding no obvious decline-rate 
relations in the JHK S bands, concluding that SNe la might well be nearly perfect standar d 
candles in the NIR at the ± 0.20 mag level or better. A study by IWood-Vasey et al.l (120081 ). 
which did not correct for optical light-curve shape, found SNe la to be excellent standard 
candles in the NIR, deriving an intrinsic dispersi on in absolute magn itude of only 0.28 
mag in the J-band and 0.15 mag in the if-band. iFolatelli et al.l (120101 ) found a marginal 
dependence of absolute NIR magnitudes on decline-rate using CSP SNe la (0.44 ± 0.14, 
0.58 ± 0.09, 0.33 ± 0.18 mag unit -1 decline-rate in YJH, respectively). Some of these 
earlier studies, however, employ inhomogeneous samples of SNe la (e.g., different telescopes, 
different reduction procedures), so subtle systematics may complicate interpretation of the 
results. 

In this article we re-examine the standardizability of SNe la in the NIR bands using a 
homogeneously obtained sample of 27 low-redshift objects observed by the CSP between 2004 
and 2007. We use these data to quantitatively examine whether a decline-rate dependence 
correction is needed in the NIR or whether SNe la are indeed perfect standard candles at 
these wavelengths. A detailed descri ption of the CSP o bserv ations, data red u ction , and 
photometry processes can be found in lHamuy et al.l ( 20061 ) and IContreras et al.l ( 2010l ). All 
of the photometric data are o n the Swope+CSP natu ral photometric system, with the final 
SN photometry published by IContreras et al.l (120101 ). with the exception of three SNe la 
(SNe 2006et, 2007af, and 2007on) wh ose photometry are presented in the second CSP data 



release article (IStritzinger et al.l 12011 



This article is organized as follows. In §2 we briefl y discuss light-cur ve template-fitting 
using SNooPy (SuperNova in Object-Oriented Python; iBurns et al.ll201ll ). and we point out 
difficulties with the template-fitting approach. In §3 we examine the standardizability of 
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SNe la in NIR bands by examining the efficacy of the Phillips method to fit the relationship 
between absolute peak Y, J, and H magnitudes and decline-rates. In §4 we investigate the 
bimodal distribution of the fast-declining, low-luminosity SNe at these wavelengths. We 
summarize our conclusions in 55. 



2. NIR Light-Curve Fitting and Morphology 

2.1. Light-Curve Fitting Using SNooPy 

From 2004-2007, the CSP obtained ~ 100 SN la light-curves in both optical and NIR 
bands. Here, we choose 27 of the best-observed objects discovered during the first three 
campaigns, all of which have optical observations that start before maximum light and 
continue for at least 30 days post-maximum. Thirteen of the SNe la have NIR observations 
that begin before maximum brightness as well, while the others begin shortly after maximum 
brightness. We denote a best-fit (BF) subsample as those 13 objects with observed maxima 
in both optical and NIR bands. This BF subsample is used to create a "training set" for the 
template light-curves following the process described here. 

SN la light-curve parameters are typically measured with respect to the peak of the light- 
curve. We prefer to directly measure the if-corrected, de-reddened time of maximum (tx ax ), 
peak magnitude (mJJl ax ), and decline-rate [Amisf- B)] based on a cub ic-spline interpolation 



of the data points themselves (see the discussion in lBurns et al.ll201ll ). However, in YJHK S , 
SNe la typically achieve peak brightness 3-4 days prior to S-band maximum, which often 
results in having objects with well-observed optical peaks but NIR observations that begin 
post-peak. Although not optimal, peak values (i.e., f x ax and m x ax ) can be estimated for such 
objects by fitting their post-peak data with a template light-curve based on observations of 
SNe la with well-observed optical and NIR peaks. 



To generate and fit our templates, we employ the SNooPy package (IBurns et al.ll201ll ). 
To create the training set using the best-observed objects, SNooPy examines the S-band 
light-curves and computes where the derivative is zero on the cubic-spline fit. This point 
is used to estimate tx ax and rrH x ax (and their uncertainties). Once f^ ax is found, SNooPy 
measures Amis^B) and then m max in all of the other filters. This training set of SNe la is 
th en used to cr e ate th e template ugriBVY J H light-curves through the procedure detailed 



by IBurns et al.l (120 111 ). Briefly, t max , m max , and Amxs(5) of the training set are placed on 



a three-dimensional surface for each band. SNooPy interpolates acros s the surface along a 



constant Arrii^B) line using a 2D variation of the Gloess algorithm ( jPersson et al.l 12004 ). 



From the template light-curves, SNooPy then measures the best-fit values of t max , m max , and 
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Amis, which is the template-derived value of the decline-rate derived from a combination of 
all filters, via x 2 minimization. For consistency, we use Am 15 (B), which is the decline-rate 
directly measured from the 5-band, as the decline-rate parameter for all of our objects in all 
filters, including those objects that were template fit. The measured Am 15 (5) parameters 
for all SNe are listed in Table HJ 



2.2. Light-Curve Morphology and the Second Maximum in the NIR 



Unlike optical light-curves, which are quite homogeneous for a fixed Am 15 (B), light- 
curves at longer wavelengths display a more diverse morphology, particularly due to the 
existence of a rise to a "second maximum" following the initial peak. Investigation in the 
i-band has shown that this second peak can exhibit a strength and morphology that varies 



significantly from SN la to SN la with identical Ami5(5) values (e.g.. iFreedman et al.ll2009 



Folatelli et al.ll2010l ; iBurns et al.ll201ll ). We confirm this with our sample in YJH (see Figs. 



[Hand [2]). Note the similarity in the B and V light-curves in Figures [T] and [2] and the difference 
in the YJH light-curves, especially around the second maximum. This creates a problem 
for template fitting NIR light-curves using a one-parameter descriptor of light-curve shape, 
especially for fast-declining SNe la [i.e., Ami 5 (B) > 1.7]. 

To test the accuracy of light-curve parameters derived from template fits, we compare 
the peak YJH magnitudes obtained from template fits with those obtained directly from 
spline fits for the BF group [excluding SNe with Ami^(B) > 1.7] in Figure |3j The weighted 
averages of the difference in peak magnitude for YJH are 0.02 mag, 0.02 mag, and 0.05 
mag, respectively, for objects with Am 15 (S) < 1.7. Figure [3] shows evidence of systematic 
differences, particularly in H, but the systematic errors are not large compared with the 
uncertainties in the final absolute magnitudes due to the errors in the host color excess and 
the derived distances. 



Folatelli et al.l (120101 ) also found evidence of systematic differences in peak magnitudes 



derived from template vs. spline fits for the NIR filters, as shown in their Figure 5, stating 
that the poor precision of the template fit in iYJH is partly due to the small sample used 
to derive t he templates and al so due to the variation in morphology surrounding the second 



maximum. 



Burns et al.l (120111 ) used a bootstrap technique to incorporate the extra dispersion 



found in the template fits caused by NIR light-curve variations. We have a ttempted to 



accou nt for this extra dispersion by adding the extrapolation errors derived by IBurns et al. 



( 1201 ll ) in quadrature with the uncertainties in the apparent magnitude obtained from the 
SNooPy fits. 
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Because the variations in the strength of the secondary maximum affect the accuracy 
of light-curve template fitting in the iYJH bands, we might expect the uncertainties in the 
peak NIR magnitudes to be a function of how many days past maximum the observations 
begin. iFolatelli et al.l (120101 ) found that if a SNe la has observations that begin within ~1 
week after the time of maximum, the random uncertainty in peak magnitude is ~0.1 mag and 
the systematic uncertainty is only ~0.03 mag; a template fit for a SNe la with photometry 
that starts later than this will be unreliable. Two thirds of the events in our sample with 
Am 15 (B) < 1.7 have photometry that starts within 5 days of the NIR maximum, and so it 
is interesting to see if the template fitting procedure gives reasonable estimates of the peak 
magnitude for these SNe la. To test this, we create a plot similar to that of Figure [3) For this 
test, we re-derive the template fit peak NIR apparent magnitudes by removing all of the data 
prior to 5 days after NIR maximum for each SN la and running the data through SNooPy 
again. Figure H] shows the resulting plots in YJH. There appear to be some systematic 
differences, but, again, they are not large when compared wih the uncertainties in the final 
absolute magnitudes. The weighted averages of the differences are 0.03 mag in Y, 0.01 mag 
in J, and 0.04 mag in H, which is about the same for the weighted averages found in Figure 
[31 suggesting that when the observations begin within 5 days of NIR maximum, SNooPy 
does an adequate job of deriving peak light-curve parameters from template fits. 

It is unclear at this point how to best handle the diversity of light-curve morphologies 
in the NIR when applying templates to SNe la whose NIR observations start more than 
5 days after NIR maxima. Introducing a second parameter, such as another peak magni- 
tude, decline-rate relation (such as a Am 15 -like parameter de fined for the YJH bandpasses) 
host-galaxy property , or spectral feature, might help (e . g.. iKasenl 120061 ; IWang et al.l 12009 



Foley fc Kasenll201ll ; ISullivan et al.ll2010l ; iBlondin et al.l 120 111 ), but identifying this second 



parameter is difficult with the limited number of well-observed events. With significant 
differences among the NIR light-curves of SNe la with similar Ara 15 (B) values, though, tem- 
plate fitting in iYJH based on the Ami^B) parameter alone may be subject to significant 
uncertainties. Here, we proceed to cautiously apply templates when needed, while using the 
BF subsample of objects, for which no template-fitting is needed, to check our results. 



3. Investigating the Homogeneity of Peak Luminosity of SNe la in the NIR 

3.1. Distance Moduli, Color Excess, and Reddening 

In this section we use our homogenous sample of light-curves to examine the precision of 
SNe la as standard candles in the NIR ban ds. We calibrate the absolute magnitudes of our 
objects using a technique first proposed by lPhillips et al.l (119991 ). which examines the corre- 
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lation between reddening-corrected, absolute peak magnitudes vs. decline-rate according to 
the model: 

Hx = m x - M x (0) - b x [Am l5 (B) - 1.1] - R X E(B - V), (1) 



where \i x is the distance modulus, m x is the peak apparent magnitude corrected for Galactic 
reddening given in band X, M x (0) is the peak absolute magnitude for Am 15 ( J B)=l.l and 
zero color excess, b x is the slope of the luminosity vs. decline-rate relation, R x is the total- 
to-selective absorption coefficient given in band X, and E(B-V) is the SN color excess due 
to host-galaxy dust. 

To calculate the distance modulus, which is used to transform rest-frame apparent 
magnitude to absolute magnitude, we employ the following approximation to the luminosity 
distance: 



4(2CMBi H , fl M , fl\) 



(1 + Z heHo ) C 1 



CMB\ ) 



(2) 



where Zheiio is the heliocentric redshift of the host-galaxy; zqmb is the redshift of the host- 
galaxy in the cosmic microwave background (CMB) rest frame; and the standard cosmolog- 
ical p arameters are Hp = 72 km s _1 Mpc -1 (IFreedman et al.ll200ll ). Qm = 0.28, and Q\ = 
0.72 (jSpergel et al.l 120071 ). Values of Zheiio and zcmb are given in Tabled) The uncertainty 
in the velocity of the host-galaxies due to peculiar velocity is assumed to be <r 2 =0.001 (300 
km s _1 ). 

For SNe that are not in the smooth Hubble flow (i.e., z < 0.01), distances derived using 
equation (2) can be inaccurate. Seven of our objects fall into this category: SNe 2005am, 
2005ke, 2006D, 2006X, 2006mr, 2007af, and 2007on. Four of these (SNe 2005ke, 2006X, 
2006mr, and 2007on) have published Cepheid or surface brightness fluctuation distances 
for their host-galaxies, which were used (see Table [2]); the other three objects were excluded 
from our sample, leaving a sample of 24 SNe la. 

For our complete sample, we employ a mix of cubic-spline (when peak is observed) 
and template (when necessary) fits to the light-curves to derive the parameters of m max 
and Ami5(B). Table [1] lists the fitting method for each SN la, and Table [3] lists the it- 
corrected spline and template-derived apparent peak magnitude in BVYJH for each object, 
it-c orrections are appl ied to convert the observed magnitudes to rest-frame magnitudes using 
the 



sec 



Hsiao et al 



Burns et all (1201 lh . 



fj2007h spectral templates. For further details on the it-correction procedure, 



All peak mag nitudes have been co rrected for Galactic reddening using the values E(B — V) ga i 
given in Tabled] (ISchlegel et al.lll998l ) and adopting Ry al =3.1 (a further discussion on this 
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point follows). In order to estimate extinction in the host-galaxies, we have used color ex- 
cesses E(B — V) obtained from the observed SNe colors as follows. For SNe with Am 15 (5) 
< 1.7, the host-galaxy color excesses in column (7) of Table [Tj are der i ved u sing the intrin- 
sic B max - V max color law at maximum light derived by iFolatelli et al.l (120101 ) . The Folatelli 
et al. relation is only valid for Ami 5 (B) < 1.7; therefore, to calculate EfB -V)host for the 



three fast-declining events, the re-derived CSP Lira law (1995) described by IFolatelli et al. 



( 120101 ) is used. We do not apply any priors to the color excess measurements; therefore, 
negative reddenings are possible and, indeed, expected due to the measurement uncertainty 
and intrinsic color dispersion in SNe la. 

Significant debate exists on the "best" value of the total-to-selective absorption coef- 
ficient for the host-galaxy, R v , for SNe la studies. By reducing the scatter in the Hub- 



ble diagram, some studies fe.g.. iTripp fc Branch 



2007 



19991: IWanget al.lbood : [Astier et al.l 12006 



Nobili k Goobar 



2008|) find that SNe la "prefer" 



Elias-Rosa et al.ll2006l ; iConley et al.l 
a lower Ry value (i?y=1.0-2.4) compared with the value obtained th r ough typical lines of 



sight in the Milky Way (R v =3.1) Oth er studies flWang et al.l 120091 : IFolatelli et al.l 12010 



Foley fc Kasenll201ll ; iMandel et al.ll2010l ) that examine the colors of SNe la as a function of 
wavelength find that the choice between a high or low Ry value depends on the reddening 
of the particular SN. It appears that SNe la that suffer from minimal— to— moderate ex- 
tinction tend to favor a value close to i?y=3.1, while those that suffer from high extinction 
favor a lower value. These studies conclude that the low derived Ry values are, in all likeli- 
hood, not physically due to changes in the actual reddening coefficient, but rather to some 
as yet-unknown p arameter affecting the intrinsic color of SN la light-curves (see, however, 
Wang et al.l 120091 ). 



The reason for these differing Ry values is not well understood, and the exact value 
depends on the method used to derive it. Given our specific goal of removing the effects of 
dust extinction from SN la light-curves in the most accurate (i.e., physically understood) 
manner possible, we choose to use the typical Galactic value of Ry =3.1. We emphasize that 
our primary purpose in this study is to examine the dependence of intrinsic (i.e., corrected 
only for extrinsic effects, such as dust extinction) peak absolute magnitude on decline-rate 
in the NIR bands - not to minimize dispersion in the Hubble diagram or to calculate cos- 
mological parameters; thus, choosing Ry =3.1 is warranted. As we shall see, since the effects 
of dust extinction in the N IR are minimal, this cho ice i s not crucial for the re sults obtained. 
Indeed, as pointed out by Krisciunas et al.l (120001 ) and IFolatelli et al.l (120101 ) . this is one of 
the significant advantages of working in the NIR. 
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3.2. Luminosity Versus Decline- Rate Relation 

In this section we investigate the correlation between peak luminosity and decline-rate 
for 5 different subsamples (see Table H]). We investigate the 5 subsamples in order to isolate 
the possible systematic effects that individual SNe la may have on the fit. Since each Rx 
(being a function only of Ry) is fixed, we solve only for M x {0) and b x using equation (1) 
in the following analyses. 



Subsample 1. Includes all 24 SNe la (i.e., the complete sample of 27 SNe la minus the 
three for which distances are not available [see §3.1]). As shown in Tableland displayed in 
Figure all three of the NIR bands show a statistically significant luminosity vs. decline-rate 
correlation, with a dependence in all three NIR bands comparable with those in the £>-band 
found by Prieto et al. (2006). However, as seen from the next subsample, the fast-declining 
and highly reddened objects have a very large effect on the derived slope. 

Subsample 2. Excludes the three fast- declining SNe la and the two most highly reddened 
objects (SNe 2005A and 2006X). The highly reddened objects were excluded in order to 
minimize the syste matic uncertainties due to the apparent peculiar nature of the reddening 
that affects them (jFolatelli et al.ll2010[ ). We choose to remove the SNe la with Ara 15 (B) > 
1.7 because this decline-rate cu toff is widely used i n other peak luminos ity vs. decline-rate 



studies (Hamuy et al. 1996a; iPhillips et al.l Il999t Kr i sciuna s 2004a,c; iPrieto et al.l 12006 



Folatelli et al.l 120101 ). and as found by iKrisciunas et al.l (120091 ) . such objects appear to obey 
a different peak luminosity vs. decline-rate relation compared with normal SNe la and 
should therefore be modeled separately. Omitting the highly reddened and fast-declining 
events creates a shallower slope on the luminosity vs. decline-rate relation compared with 
subsample 1 for all three NIR bands (Fig. [6]), with the rms dispersion decreasing significantly 
in all bands as well. There is evidence for a weak correlation between peak luminosity and 
decline-rate in Y, with the significance of the measured slope being 1.5a, and a strong 
correlation in both J and H at the 3a significance level. 

Subsample 3. Uses SNe la, with first observations starting within 5 days after NIR peak 
brightness. This subsample also omits the highly reddened and fast-declining events. Like 
subsample 2, the Y-band shows a marginal (1.7cr) dependence of absolute peak magnitude 
on decline-rate, while both the J- and //-bands show a stronger (2.8a) correlation between 
peak absolute magnitude and decline-rate (see Tableland Fig. [7]). As shown in Figure [7] 
and Table |5l the rms scatter decreases significantly in YJH, compared with subsample 2. 



Subsample 4- Uses the BF subsample, which ensures that purely empirical values for the 
peak magnitude and Am 15 (5) are used, and avoids the unknown systematics that plague 
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the template fitting process. This set produces the largest peak luminosity vs. decline-rate 
dependence of any set examined (see Fig. [8] and Table [5]), but as shown by subsample 5, this 
is likely due to the inclusion of the highly reddened SN 2006X and the three fast-declining 
events. 

Subsample 5. Same as subsample 4, but excludes the highly reddened SN 2006X and 
the fast-declining [Am^S) > 1.7] SNe la. This subsample has 7 objects with a maximum 
Aro 15 (5) of 1.39. For this restricted sample, there appears to be only little dependence 
of peak luminosity on decline-rate within the uncertainties for the Y- and J-bands. The 
calculated slopes are 0.14 ± 0.47 and 0.31 ± 0.26 for Y and J, respectively. In H, there is a 
marginal (2a) detection of a correlation between peak luminosity and decline-rate (see Fig. 
IH]). These results are consistent with the luminosity vs. decline-rate relations derived from 
subsamples 2 and 3. The scatter on the corrected magnitudes range between 0.05-0.14 mag, 
with the smallest dispersion in H and the largest in Y. 



3.3. Comparison with Previous Studies 



To test whether corrections for decline-rate substantially improve the precision with 
which SNe la can be used as standardizable candles in the NIR, one can examine the dis- 
persion in derived absolute magnitude before and after such corrections are applied. From 
examination of Table |6j which displays both the uncorrected (bx = 0) and corrected rms 
values given an Ry = 3.1 for the 5 subsamples used throughout the article, we find that a 
significant rms decrease occurs for the J/-band, followed by the J-band. The degree of confi- 
dence for a non-zero slope in each bandpass [i.e., that a relationship exists between Am 15 (B) 
and absolute brightness] is indicated in column (6) of Table 6. Using the t-test statistic, and 
based on the degrees of freedom for the fit, the confidence level for a non-zero slope was 
calculated. For all but subsample 5, we can exclude a zero slope at the 93% confidence level 
or better. This decrease in the rms and the high confidence levels confirm our earlier findings 
that correcting observed SNe la values in the J- and J/-bands may substantially improve 
their utility as cosmological distance indicators. 

To compare with previous studies, we focus on the uncorrected dispersion in Table [6] 
for subsample 2. We note that the intrinsic (i.e., uncorrected) dispersions found for sub- 
sample 2 (0.16 mag and . 17 m ag for J and H, respectively) compare well with the rms 
values of iKrisciunas et al.l (120091 ) for SNe la whose N IR light-curves pea k befo re the epoch 
of B max (0.16 mag in J and 0.15 mag in H) and of IWood-Vasey et al.l ( 120081 ) in H (0.15 
mag). Our J-band rms y alue is much less than the J-band rms value of 0.33 mag found by 
Wood-Vasev et al.l (l2008f ). 
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3.4. The Effect of Ry on Derived Correlations 

Given that the value for Ry is a hot topic of debate in SNe la studies, we examine the 
effect that lowering Ry (from R v =3.1 to Ry =1.7) has on the slope of the peak absolute 
YJH magnitudes vs. decline-rate relation. The same five subsamples described in §3.2 are 
used for this investigation. The measured slopes are given in the first two rows for each filter 
and subsample in Table El 

The derived slopes for the two Ry values in subsample 1 are within the measured 
uncertainties for each NIR band. Subsample 2 produces the smallest variance between 
slopes, given different values of Ry. Subsamples 3 and 4 show a larger deviation between 
slopes in all three NIR bands, but all slopes are within measured uncertainties. The derived 
slopes for subsample 5 show a smaller variance than subsamples 1, 3, and 4. For all five 
subsamples, the slopes are rather insensitive to the exact value of Ry. 



3.5. Hubble Diagram 

Using the distance modulus derived from equation (1) and the CMB redshifts taken 
from NASA/IPAC Extragalactic Database (NED), we construct a Hubble diagram, which is 
plotted in Figure (TDJ We use SNe and fit parameters from subsample 2 to derive the distance 
modulus for each SN in YJH and 5-bands \M B (0) = -19.222(005), b B = 1.042(213)]. The 
distance moduli for the three NIR bands are averaged for each object. The resulting Hubble 
diagram for the combined NIR bands and S-band is shown in the top panel of Figure [TUl The 
solid line represents the adopted concordance model of equation (2). The residuals between 
the corrected distance modulus and the standard cosmology are shown in the bottom panels 
of Figure [10J We find a combined rms scatter of 0.13 mag, or 6% in distance, for the NIR 
bands and a rms scatter of 0.22 mag, or 11% in distance, for the S-band. The decrease in 
scatter from optical bands to NIR bands on the Hubble diagram provides strong evidence 
that SNe la in NIR bands are excellent standard candles. 

Note that the dispersion of 0.13 mag obtained in the N IR includes the effects of peculiar 



velocities. This was demonstrated by lFolatelli et al.l (120101 ). who pointed out that the resid- 



uals in the Hubble diagram in separate filters are highly correlated. Peculiar velocities also 
explain why the dispersion in Figure [TUl decreases with redshift. Hence, the true precision of 
SNe la as distance indicators in the NIR is almost certainly better than 6%. 
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4. Bimodal Distribution of Absolute Magnitudes for Fast-Declining SNe la 



Figures [5] an d [8] indicate a range of peak luminosities for the fast-declining SNe la. 
Krisciunas et all (120091 ) suggest that, instead of a range of peak brightnesses, there is a 
bimodal distribution of peak absolute NIR magnitudes for the fast-declining objects; those 
with NIR maximum occurring after the time of 5-band maximum (fg ax ) are sub luminous 
and those with NIR maximum occurring prior to f^ ax have standard brightnesses compared 
with "normal" SNe la. Our data are consistent with this bimodal distribution hypothesis. 
Our sample has three fast-declining SNe la: SN 2005ke (Amxs(S) = 1.75), SN 2006mr 
(Am 15 ( J B)=1.81), and SN 2007on (Am 15 (S)=1.89). As shown in Table El SNe 2005ke and 
2006mr have an average NIR maxima that occurs ~2-4 days after t£ ax and are, on average, 
0.89 mag fainter in My, 0.98 mag fainter in Mj, and 0.63 mag fainter in M h than the average 
peak luminosity of the slow to midrange-decliners, respectively, whereas SN 2007on has NIR 
maximum that occurs ~2 days prior to f§ ax and has standard peak absolute brightness 
in all three NIR b ands. It is not s urpris ing that this same bimodality in the absolute 
magnitude found by lKrisciunas et all (120091 ) is also found i n our sample seeing, as o ur sample 
contains some of t he same objects that we re used in the iKrisciunas et all (120091 ) study. It 
was also found by IKrisciunas et all (120091 ) that fast-decliners with standard brightnesses, 
such as SN 2007on, show a distinct secondary maxima in the NIR, much like slow- declining 
objects, whereas the subluminous fast-decliners have monotonically declining light-curves in 
all bands. Our three fast-declining objects show this trend (see Figure [2]). 

We re-examine the peak absolute magnitude vs. decline-rate relationship, excluding 
the two subluminous fast-decliners and including the standard fast-decliner. Subsamples 2 
and 5 described in §3.2 are used for this analysis, SNe 2005ke and 2006mr are omitted and 
SN 2007on is included for each subsample. The fit parameters derived are given in Table [SJ 

Including SN 2007on and excluding the two subluminous fast-decliners creates no signif- 
icant change in slopes for the three NIR bands in subsamples 2 and 5, with the significance 
on the measured slopes remaining the same. The only case where the fiducial absolute peak 
magnitude significantly changes when SN 2007on is included is for the J-band in subsample 
2. The absolute peak magnitudes fo r the other band s in su bsamples 2 and 5 change very 
little. Our findings here confirm the IKrisciunas et al.l (120091 ) results that SNe la with NIR 
maximum occurring prior to i?-band maximum, such as SN 2007on, have absolute peak 
magnitudes that are consistent with the absolute peak magnitudes of slow- to mid-decliners, 
allowing these SNe la to be used for the peak absolute magnitude vs. decline-rate relationship 
analysis. 
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5. Discussion and Conclusion 

We analyze 27 of the best-observed SNe la from the first 4 years of the CSP. Spline 
curves were fit to the 13 objects with observed maxima in YJH. A "training set" was 
created from these 13 SNe la to construct light-curve templates for the other 14 objects 
with no observed NIR maxima. Light-curve properties obtained from these fits were used to 
calculate peak absolute magnitude, and correlations were sought with £>-band decline-rates. 

We discuss difficulties encountered when fitting templates to the NIR light-curves, due 
to the variations in the strength of the secondary maximum at a given decline-rate. The 
problems are especially pronounced for fast-decliners; however, if the observations begin 
within 5 days after NIR maximum, SNooPy provides peak magnitudes with an accuracy 
of ~0.04 mag. Obtaining additional SNe la light-curves with observed NIR maxima would 
be ideal, since splines can be directly fit, but this poses an observational challenge, since 
NIR maximum typically occurs several days earlier than it does in optical bands. Finding a 
second light-curve parameter [i.e., other than Am 15 (B)] could certainly be beneficial when 
fitting NIR light-curve templates to data; our limited sample size, however, makes such 
an investigation problematic. Future larger samples of SNe la observed in the NIR could 
possibly provide evidence for this second parameter and increase the accuracy of fitting 
templates to NIR light-curves, regardless of when the observations begin. 

From our original sample of 27 objects, we examine whether a "Phillips relation" exists 
for the observed data that improves the standardizability of SNe la in the NIR for five 
different subsamples. Ry was kept fixed at 3.1, although the results are shown to not change 
appreciably if lower values are used. Consistent with earlier findings, we find that any peak 
luminosity vs. decline-rate dependence does indeed decrease in the NIR compared with the 
optical, confirming that SNe la are better "standard candles" at longer wavelengths. 

When fast-decliners and highly reddened events are removed from the sample (e.g., 
subsample 2), a weaker dependence of absolute peak magnitude on decline-rate is found, 
with the significance of the measured slope being 1.5a for the F-band and 3a for the J- and 
if-bands. The fact that the rms dispersion decreases significantly when the fast-declining 
SNe la are removed confirms the finding by Krisciunas et al. (2009) that a simple linear 
relationship cannot be used to fit both the slow- to midrange-decliners and the fast-decliners. 
For subsample 3, which uses SNe la whose first observations begin within 5 days after NIR 
maxima, we find a weak (1.7a) dependence of Y peak luminosity on decline-rate and a 
stronger (2.8a) dependence of J- and H-peak luminosities on decline-rate. The dispersion 
decreases significantly for subsample 3, compared with subsample 2. Subsample 4, which 
uses the BF subsample, shows the largest correlation between peak absolute magnitude and 
decline-rate, but the fast-declining SNe la strongly affect the slope. For subsample 5, which 
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uses the BF subsample excluding the fast-decliners and the highly reddened object, there 
is no significant dependence of peak Y luminosity on decline-rate. There appears to be a 
weak (1.2a) correlation between decline-rate and the J-band absolute peak magnitude in 
subsample 5 and still a marginal (2.0a) correlation between decline-rate and absolute peak 
magnitude in H. After removal of the highly reddened and fast-declining SNe la, we find 
evidence at the 2. 0-3. Oct level for a luminosity vs. decline-rate relations in J and H, but 
minimal evidence for such a relation in Y. 

The lowest dispersions in the uncorrected absolute magnitudes are obtained when the 
highly reddened and fast-declining SNe la are removed from the sample. Comparing the 
uncorrected and corrected values, we are at least 93% confident that the slope is nonzero 
for subsamples 1-4. It could be beneficial to correct the peak luminosity in both the J- 
and .ff-bands in order for SNe la to be utilized as precise standard candles in the NIR. The 
dispersions obtained are smaller than those found in previous NIR studies. 

As larger samples of SNe la in the NIR are gathered, and especially as more objects 
are observed pre-maximum in the NIR, an interesting study could be to investigate whether 
using a Am 15 -like relation, defined in the NIR bands (rather than in the £>-band) would 
further decrease the dependence of absolute peak magnitude on decline-rate. Our sample 
presented here does not include enough objects that have observed NIR maxima to test this 
idea, but future SNe la studies will certainly have more observations that cover the NIR 
peak with greater frequency. 

Evidence for a bimodal distribution was found for the NIR absolute peak magnitudes of 
our fast-declining SNe la. The SN from our sample that peaked before the time of £?-band 
maximum (SN 2007on) had a "normal" peak luminosity, whereas those that peaked after 
the time of £>-band maximum (SNe 2005ke and 2006mr) had subluminous absolute peak 
magnitudes compared with the rest of the sample. It will be interesting to see if future 
SNe la studies find any slow- to midrange-decliners that possess this bimodal distribution. 
It would also be beneficial to obtain larger samples of the fast-declining SNe la at NIR, in 
addition to optical passbands, to further our understanding of these "peculiar" SNe la. 

We also present a Hubble diagram for the NIR bands and the 5-band, using the pa- 
rameters derived from subsample 2. The rms scatter greatly decreases going from optical 
to NIR bands (0.22 mag in B to 0.13 mag for combined NIR bands), again demonstrating 
that SNe la observed in NIR bands are excellent standard candles. The Hubble diagram and 
Hubble residuals in Figure [TU| along with the weaker dependence of absolute peak magnitude 
on decline-rate found for the NIR bands compared with the optical bands, clearly show the 
advantages of working in the NIR for distance determinations and cosmological studies. Fu- 
ture high-z SNe la cosmological studies may therefore benefit from observations taken in the 
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rest-frame NIR bands, since SNe la appear to be better standard candles in these passbands 
than in the optical (i.e., they have a weaker dependence of absolute NIR peak magnitude 
on decline-rate, systematic uncertainties are minimized due to the reduced effects of dust 
extinction, and NIR passbands are relatively insensitive to the exact value of Ry used). The 
y-band, in particular, appears to hold significant promise for such future studies, since it 
offers an optimal combination of signal-to-noise ratio, low sensitivity to dust extinction, and 
an essentially negligible dependence of absolute peak magnitude on decline-rate. 
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Table 1. List of Supernovae and Properties 



G1\T 


Zhelio 


ZCMB 




h{B - V) ga i 


- V J host 


• 

Fit method 


Subsample 


I 1 ) 


\ 2 ) 


(3) 


( 4 ) 




(6) 


(i\ 
{') 




(8) 


2004ef 


0.0310 


0.029 


1.345(017) 


0.056(006) 


0.097(008) 


Template 


1, 


2,3 


2004eo 


0.0157 


0.0147 


1.370(055) 


0.108(011) 


0.014(018) 


Spline 


1, 


2, 3, 4, 5 


2004ey 


0.0158 


0.0146 


0.935(014) 


0.139(014) 


-0.052(005) 


Template 


1, 


2, 3 


2004gs 


0.0267 


0.0275 


1.502(016) 


0.031(003) 


0.157(005) 


Template 


1, 


2 


2004gu 


0.0459 


0.0469 


0.900(045) 


0.026(003) 


0.183(009) 


Template 


1, 


2,3 


2005A 


0.0191 


0.0184 


1.224(102) 


0.030(003) 


0.996(016) 


Template 


1 




2005M 


0.0220 


0.0230 


0.859(011) 


0.031(003) 


0.052(011) 


Spline 


1, 


2, 3, 4, 5 


2005ag 


0.0794 


0.0801 


0.905(046) 


0.041(004) 


0.047(007) 


Template 


1, 


2,3 


2005al 


0.0124 


0.0133 


1.146(015) 


0.055(005) 


-0.094(005) 


Template 


1, 


2,3 


2005am 


0.0079 


0.0090 


1.509(021) 


0.054(005) 


-0.030(010) 


Template 






2005el 


0.0149 


0.0149 


1.345(020) 


0.114(011) 


-0.043(016) 


Spline 


1, 


2, 3, 4, 5 


2005eq 


0.0290 


0.0284 


0.756(022) 


0.074(007) 


0.090(006) 


Template 


1, 


2,3 


2005hc 


0.0459 


0.0450 


0.894(014) 


0.029(003) 


0.048(006) 


Template 


1, 


2,3 


2005iq 


0.0340 


0.0330 


1.258(013) 


0.022(002) 


-0.040(006) 


Template 


1, 


2,3 


2005kc 


0.0151 


0.0139 


1.163(029) 


0.132(013) 


0.173(006) 


Spline 


1, 


2, 3, 4, 5 


2005ke 


0.0049 


0.0045 


1.748(015) 


0.027(003) 


0.036(005) 


Spline 


1, 


, , 4 


2005ki 


0.0192 


0.0204 


1.279(038) 


0.032(003) 


-0.018(012) 


Spline 


1, 


2, 3, 4, 5 


2006D 


0.0085 


0.0097 


1.388(012) 


0.046(005) 


0.071(010) 


Spline 






2006X 


0.0052 


0.0063 


1.107(019) 


0.026(003) 


1.202(013) 


Spline 


1 




2006ax 


0.0167 


0.0180 


0.990(016) 


0.050(005) 


-0.033(006) 


Spline 


1, 


2, 3, 4, 5 


2006bh 


0.0109 


0.0105 


1.426(011) 


0.026(003) 


-0.038(007) 


Template 


1, 


2, 3 


2006eq 


0.0495 


0.0484 


1.660(060) 


0.048(005) 


0.191(033) 


Template 


1, 


2 


2006et 


0.0222 


0.0212 


0.903(012) 


0.019(002) 


0.222(008) 


Spline 


1, 


2, 3, 4, 5 


2006gt 


0.0448 


0.0437 


1.663(037) 


0.037(004) 


0.165(019) 


Template 


1, 


2 


2006mr 


0.0059 


0.0055 


1.805(015) 


0.021(002) 


-0.025(013) 


Spline 


1, 


, ,4 


2007af 


0.0055 


0.0063 


1.203(036) 


0.039(004) 


0.177(019) 


Spline 






2007on 


0.0065 


0.0062 


1.893(018) 


0.011(001) 


0.022(010) 


Spline 


1, 


, ,4 
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Note. - Col. (1) SN name; Col. (2) Heliocentric redshift from the NASA/IPAC Extragalactic 
Database (NED); Col. (3) Redshift in the frame of the 3K CMB (NED); Col. (4) Observed 
Ami5 (.B); Col. (5) Galactic reddening; Col. (6) Host-galaxy reddening; Col. (7) Fitting method 
(spline or template); Col. (8) Subsample(s) SN belongs to (see text and Table 4) If blank, SN is 
excluded from all subsamples. 
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Table 2. Derived Distance Moduli 



SN 


Host-galaxy 




(mag 


;) 


Distance method 


Reference 


(1) 


(2) 




(3) 




(4) 


(5) 


2004ef 


UGC 12158 




35.56 


± 


0.07 


LD 




2004eo 


NGC 6928 




34.00 


1 

± 


0.15 


LD 




2004ey 


UGC 11816 




33.99 


± 


0.15 


LD 




2004gs 


MCG +03-22- 


020 


35.31 


± 


0.08 


LD 




2004gu 


FGC 175A 




36.50 


1 


0.05 


LD 




2005A 


NGC 958 




34.47 


± 


0.12 


LD 




2005M 


NGC 2930 




34.94 


± 


0.10 


LD 




2005ag 






37.72 


i 

± 


0.03 


LD 




2005al 


NGC 5304 




33.71 


± 


0.16 


LD 




2005am 


NGC 2811 




32.85 


l 


0.24 


LD 




2005el 


NGC 1819 




33.99 


1 

± 


0.14 


LD 




2005eq 


MCG -01-09- 


006 


35.43 


± 


0.08 


LD 




2005hc 


MCG +00-06- 


003 


36.47 


± 


0.05 


LD 




2005iq 


MCG -03-01- 


008 


35.78 


± 


0.07 


LD 




2005kc 


NGC 7311 




33.88 


± 


0.16 


LD 




2005ke 


NGC 1371 




31.84 


± 


0.08 


SBF 


Tonrv et al. ( 20011* 


2005ki 


NGC 3332 




34.64 


± 


0.11 


LD 






MCG —01 — ^H- 


OH 


33.00 


± 








2006X 


NGC 4321 




30.91 


± 


0.14 


Cepheids 


Freedman et al. (2001) 


2006ax 


NGC 3663 




34.36 


± 


0.12 


LD 




2006bh 


NGC 7329 




33.23 


± 


0.21 


LD 




2006eq 






36.64 


± 


0.05 


LD 




2006et 


NGC 232 




34.80 


± 


0.10 


LD 




2006gt 


2MASX J00561810- 


0137327 


36.41 


± 


0.05 


LD 




2006mr 


NGC 1316 




31.15 


± 


0.23 


SBF 


Aihar et al. f 20011 


2007af 


NGC 5584 




32.10 


± 


0.34 


LD 




2007on 


NGC 1404 




31.45 


± 


0.19 


SBF 


Jensen et al. (2003) 
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Note. - Col. (1) SN name; Col. (2) Host-galaxy; Col. (3) Distance modulus from fits of §4; Col. 
(4) Distance method used (LD - luminosity distance, SBF - surface brightness fluctuation); Col. (5) 
Reference used for SBF and Cepheid distances. 



*NGC 1371 is a member of the Eridanus group/cluster. iTonry et al.l (120011 ) give a distance modulus 
of 32. 00 ± 0.08 mag for this group based o n SBF measurem e nts o f seven members. Subtracting 0.16 
mag (I Jensen et al.l 120031 ) to put this on the iFreedman et al.l (j200lf ) scale gives 31.84 ± 0.08 mag. 
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Table 3. Supernova Light-Curve Properties 



Filter 


tmax (JD -2,453,000) 


m max (spline) 


m max (template) 




(days) 


(mag) 


(mag) 


(1) 


(2) 


(3) 


(4) 



SN 2004ef 



B 


264 


40 


± 





05 




17.077 


± 





007 


V 


265 


88 


± 





05 




16.910 


± 





004 


Y 


261 


47 


± 





05 




16.999 


± 





008 


J 


261 


44 


± 





05 




17.020 


± 





114 


H 


261 


04 


± 





05 


SN 2004eo 


17.409 


± 





027 


B 


278 


69 


± 





28 


15.497 ± 0.013 


15.552 


± 





Oil 


V 


280 


68 


± 





54 


15.359 ± 0.012 


15.346 


± 





006 


Y 


275 


85 


± 





91 


15.929 ± 0.102 


15.863 


± 





027 


J 


274 


85 


± 





19 


15.547 ± 0.010 


15.515 


± 





025 


H 


274 


18 


± 





13 


15.814 ± 0.010 


15.742 


± 





024 


B 


304 


87 


± 





04 


SN 2004ey 


15.275 


± 





004 


V 


306 


20 


± 





04 




15.224 


± 





003 


Y 


300 


61 


± 





04 




15.747 


± 





033 


J 


301 


40 


± 





04 




15.540 


± 





021 


H 


300 


88 


± 





04 




15.757 


± 





018 


B 


356 


23 


± 





07 


SN 2004gs 


17.277 


± 





003 


V 


357 


92 


± 





07 




17.057 


± 





004 


Y 


353 


77 


± 





07 




17.241 


± 





022 


J 


353 


81 


± 





07 




17.138 


± 





032 


H 


353 


52 


± 





07 




17.236 


± 





032 


B 


362 


56 


± 





20 


SN 2004gu 


17.570 


± 





007 


V 


363 


91 


± 





20 




17.401 


± 





005 


Y 


357 


98 


± 





20 




17.991 


± 





046 
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Table 3 — Continued 



r liter 


t max (JD -2,453,000) 


m max (spline) 


m max (template) 




(days) 


(mag) 




(mag) 


(1) 


(2) 




(3) 




(4) 


j 


359.01 ± 


0.20 




17. 


.934 


± 0.093 


ff 


358.28 ± 


0.20 


SN 2005 A 


18. 


403 


± 0.050 


R 


378.86 ± 


0.24 




18. 


.239 


± 0.013 




380.18 ± 


0.24 




17. 


.214 


± 0.010 


Y 


374.59 ± 


0.24 




16 


,306 


± 0.022 


j 


375.37 ± 


0.24 




16 


429 


± 0.023 


ff 


374.85 ± 


0.24 


SN 2005M 


16 


,302 


± 0.022 


ft 


406.02 ± 


0.20 




16 


,048 


± 0.003 


V 


407.49 ± 


0.12 


16.003 ± 0.010 


15. 


.998 


± 0.002 


Y 


402.52 ± 


0.81 


1 6 627 4- 036 


16 


.587 


± 0.006 


7 


402.59 ± 


0.94 


1 6 459 + 046 


16 


.463 


± 0.013 


ff 


401.38 ± 


0.06 


16 4^8 4- 074 
SN 2005ae- 


16 


6^2 


+ 01 


B 


414.79 ± 


0.14 




18. 


.607 


± 0.005 


V 


416.22 ± 


0.14 




18. 


.558 


± 0.005 


Y 


410.17 ± 


0.14 




19. 


,325 


± 0.032 


J 


411.10 ± 


0.14 




19. 


,244 


± 0.048 


H 


410.53 ± 


0.14 


SN 2005al 


19. 


,044 


± 0.063 


B 


429.86 ± 


0.13 




15 


,045 


± 0.003 


V 


431.11 ± 


0.13 




15. 


.094 


± 0.004 


Y 


426.42 ± 


0.13 




15. 


.576 


± 0.035 


J 


426.55 ± 


0.13 




15. 


,355 


± 0.048 


H 


426.04 ± 


0.13 


SN 2005am 


15. 


.706 


± 0.021 


B 


435.90 ± 


0.13 




13 


,870 


± 0.008 
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Table 3 — Continued 



b liter 


t max (JD -2,453,000) 
(days) 


m max (spline) 
(mag) 


m max (template) 
(mag) 


(1) 


(2) 




(3) 


(4) 


V 

V 


437.47 ± 


0.13 




13.750 


± 0.006 


Y 


433.30 ± 


0.13 




13.949 


± 0.032 


j 


433.29 ± 


0.13 




13.746 


± 0.045 


ff 

±1 


432.97 ± 


0.13 


SN 2005p1 


14.103 


± 0.022 


ft 


647.30 ± 


0.36 


1 5 278 + 011 


15.307 


± 0.012 


V 


647.86 ± 


0.24 


15.194 ± 0.012 


15.261 


± 0.005 


Y 

1 


642.66 ± 


0.15 


15 559 + 014 


15.676 


± 0.012 


j 


643.15 ± 


0.15 


1 5 508 4- 008 

± *J . tJ W(J _1_ U.uu(j 


15.574 


± 0.015 


ff 


642.33 ± 


0.12 


15 652 + 01 


15.697 


± 0.012 




654.60 ± 


0.10 




16.565 


± 0.004 


V 

V 


655.93 ± 


0.10 




16.458 


± 0.004 


Y 


650.09 ± 


0.10 




1 7 062 


4- 020 


j 


651.11 ± 


0.10 




16.868 


± 0.019 




650.39 ± 


0.10 




17.122 


± 0.022 




667.71 ± 


0.07 




17.431 


± 0.004 


V 


669.10 ± 


0.07 




17.395 


± 0.004 


Y 


663.18 ± 


0.07 




18.049 


± 0.025 


J 


664.14 ± 


0.07 




17.922 


± 0.026 


H 


663.55 ± 


0.07 




18.003 


± 0.051 


B 


688.03 ± 


0.05 


SN 2005iq 


16.887 


± 0.004 


V 


689.29 ± 


0.05 




16.902 


± 0.004 


Y 


684.42 ± 


0.05 




17.442 


± 0.021 


J 


684.61 ± 


0.05 




17.365 


± 0.016 


H 


684.08 ± 


0.05 




17.513 


± 0.025 
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Table 3- 


—Continued 




Filter 


t max (JD -2,453,000) 


m max (spline) 


m max (template) 




(days) 


(mag) 


(mag) 


(1) 


(2) 


(3) 


(4) 









SN 2005kc 












B 


698.29 ± 


0.22 


16.037 ± 0.004 


16 


.063 


± 


0. 


.012 


V 


699.19 ± 


0.08 


15.744 ± 0.004 


15. 


.745 


± 


0. 


006 


Y 


693.97 ± 


0.17 


15.658 ± 0.014 


15. 


,626 


± 


0. 


.014 


J 


694.69 ± 


0.73 


15.479 ± 0.010 


15. 


.455 


± 


0. 


016 


H 


694.10 ± 


0.55 


15.684 ± 0.013 


16 


,644 


± 


0. 


.010 








SN 2005ke 












B 


698.81 ± 


0.07 


14.882 ± 0.006 


14 


,890 


± 


0. 


013 


V 


701.32 ± 


0.24 


14.189 ± 0.005 


14 


,202 


± 


0. 


.018 


Y 


701.26 ± 


0.70 


14.055 ± 0.018 


14. 


,078 


± 


0. 


.008 


J 


700.15 ± 


0.10 


14.048 ± 0.006 


14. 


.016 


± 


0. 


.010 


H 


702.07 ± 


0.13 


13.932 ± 0.008 


13 


.928 


± 


0. 


004 








SN 2005ki 












B 


706.13 ± 


1.28 


15.689 ± 0.011 


15. 


.690 


± 


0. 


009 


V 


706.40 ± 


0.19 


15.670 ± 0.005 


15 


,698 


± 


0. 


.007 


Y 


701.62 ± 


0.06 


16.233 ± 0.011 


16 


.278 


± 


0. 


.012 


J 


702.73 ± 


0.26 


16.129 ± 0.012 


16 


,186 


± 


0. 


020 


H 


701.48 ± 


0.75 


16.265 ± 0.028 


16 


,321 


± 


0, 


.024 








SN 2006D 












B 


757.75 ± 


0.06 


14.319 ± 0.007 


14 


,339 


± 


0. 


.008 


V 


759.13 ± 


0.19 


14.183 ± 0.007 


14. 


.212 


± 


0. 


004 


Y 


754.85 ± 


0.05 


14.754 ± 0.117 


14. 


.630 


± 


0. 


013 


J 


755.05 ± 0.067 


14.379 ± 0.008 


14. 


.460 


± 


0. 


010 


H 


754.90 ± 0.055 


14.500 ± 0.007 


14. 


.590 


± 


0. 


.008 








SN 2006X 












B 


786.23 ± 0.14 


15.324 ± 0.007 


15 


.383 


± 


0. 


021 


V 


789.18 ± 0.13 


14.111 ± 0.011 


14 


,111 


± 


0. 


.010 


Y 


781.50 ± 0.11 


13.144 ± 0.007 


13 


,175 


± 


0. 


.017 
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Table 3 — Continued 



b liter 


t max (JD -2,453,000) 


m max (spline) 


m max (template) 




(days) 




(mag) 


(mag) 


(1) 


(2) 






(3) 


(4) 


j 


782.74 ± 


0. 


12 


12 866 4- 006 

±4j . (JVJVJ _l_ W . WWVJ 


12.901 ± 0.012 


ff 


781.76 ± 


0. 


12 


12 QQfi + D D04 


12.999 ±0.010 
















R 


827.59 ± 


0. 


13 


1 Fi 201 + 00 s ) 


15.213 


± 0.007 




828.27 ± 


0. 


.36 


15.211 ± 0.002 


15.228 


± 0.005 


Y 


823.47 ± 


0. 


.57 


15 835 ± 012 


15.932 


± 0.014 


j 


823.72 ± 


0. 


.55 


15 740 + 01 1 


15.754 


± 0.012 


ff 


824.24 ± 


0. 


.40 


i r. Q22 + 021 


16.046 


± 0.010 


ft 


833.65 ± 


0. 


.04 




14.488 


± 0.006 


V 


835.00 ± 


0. 


.04 




14.477 


± 0.004 


Y 


830.57 ± 


0. 


.04 




14.978 


± 0.009 


7 


830.57 ± 


0. 


.04 




14.930 


± 0.018 


ff 


830.13 ± 


0. 


.04 




15.043 


± 0.014 


B 


976.72 ± 


0. 


.43 




18.726 


± 0.027 


V 


978.73 ± 


0. 


.43 




18.438 


± 0.019 


Y 


975.28 ± 


0. 


.43 




18.236 


± 0.018 


J 


975.60 ± 


0. 


.43 




18.247 


± 0.114 


H 


975.55 ± 


0. 


.43 


SN 2006et 


18.374 


± 0.065 


B 


993.55 ± 


0. 


.37 


16.038 ± 0.004 


16.060 


± 0.007 


V 


995.75 ± 


0. 


.21 


15.859 ± 0.008 


15.865 


± 0.004 


Y 


989.27 ± 


1. 


.07 


16.344 ± 0.016 


16.333 


± 0.013 


J 


990.48 ± 


0. 


53 


16.055 ± 0.063 


16.113 


± 0.012 


H 


989.34 ± 


0. 


.27 


16.323 ± 0.031 
SN 2006gt 


16.401 


± 0.019 



B 



1003.15 ± 0.10 



18.383 ± 0.016 
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Table 3 — Continued 



Filter 


t max (JD -2,453,000) 


m max (spline) 


m max (template) 




(days) 






(mag) 




(mag) 




(1) 


(2) 






(3) 




(4) 




V 


1005 1Q 4- 


n 


1 




1 8 


129 


± 0. 


01 1 


Y 


1001 05 + 


n 


10 




18 


326 


± 0. 


015 


J 


1 DD9 98 4- 


n 

u . 


1 n 




1 8 


,244 


± 0. 


034 


H 


1 DD9 9Q 4- 


n 
u 


1 n 




1 8 


.265 


± 0. 


0^5 










SN 2006mr 










B 


1 050 66 + 

XWtJW.VJVJ _l_ 


n 


08 


15.387 ± 0.009 


1 5 


,555 


± 0. 


01 9 


V 


1059 5^ 4- 

IUiJ^i. JO J_ 


n 
u 


(14 


14.640 ± 0.009 


14 


.780 


± 0. 


ni 6 


Y 


1 056 68 4- 


n 

u . 




14.019 ± 0.006 


14 


.075 


± 0. 


DD8 


J 




n 
u 


(19 


14.007 ± 0.004 


14 


.107 


± 0. 


ni 4 


H 


1055.30 ± 


0. 


.03 


13.850 ± 0.005 


13 


.903 


± 0. 


009 










SN 2007af 










D 
D 


1174.78 ± 


0. 


,20 


iq a An -i- n nn 

10.44U It U.U11 


13 


,450 


± 0. 


.012 


V 


1176.16 ± 





17 


13.228 ± 0.016 


13 


,235 


± 0. 


.007 


Y 


1170.57 ± 


0. 


.08 


13.582 ± 0.009 


13 


,584 


± 0. 


.007 


J 


1171.12 ± 





,04 


13.469 ± 0.002 


13 


,460 


± 0. 


.008 


H 


1170.21 ± 


0. 


.15 


13.623 ± 0.004 


13 


,625 


± 0. 


.005 










SN 2007on 










B 


1420.50 ± 


0. 


.08 


13.085 ± 0.007 


13 


.187 


± 0. 


.015 


V 


1421.78 ± 


0. 


.15 


12.973 ± 0.007 


13 


.024 


± 0. 


009 


Y 


1418.04 ± 





.04 


13.254 ± 0.005 


13 


.388 


± 0. 


.010 


J 


1418.68 ± 


0. 


.32 


13.093 ± 0.063 


13 


.239 


± 0. 


016 


H 


1418.36 ± 





.08 


13.253 ± 0.006 


13 


.285 


± 0. 


.007 



Note. - Col. (1) Filter; Col. (2) t max ; Col. (3) /("-corrected m 
spline; Col. (4) ZT-corrected m max template. 



Table 4. Description of Subsamples 



Subsample Description N$Ne 



1 All 24 SNe 24 

2 Exc. fast-declining (Am i5 (B) > 1.7) and highly reddened (E(B — V) > 0.9) 19 

3 SNe whose first observation begins < 5 days after NIR maximum 

exc. fast-declining and highly reddened events 13 

4 Best fit 11 

5 Best fit exc. fast-declining and highly reddened 7 
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Table 5. Fit Parameters of Peak NIR Magnitude Versus Am 15 (5) 



Filter X M x (0) b x RMS N SN e Subsample 

(mag) (mag decline-rate -1 ) (mag) 



(1) 


(2) 






(3) 




(4) 


(5) 


(6) 


Y 


-18.535 


± 


0.032 


0.611 


± 


0.263 


0.382 


24 


1 


J 


-18.588 


± 


0.028 


0.691 


± 


0.230 


0.336 


24 


1 


H 


-18.432 


± 


0.017 


0.581 


± 


0.153 


0.235 


24 


1 


Y 


-18.494 


± 


0.007 


0.294 


± 


0.187 


0.184 


19 


2 


J 


-18.574 


± 


0.003 


0.363 


± 


0.119 


0.141 


19 


2 


H 


-18.415 


± 


0.006 


0.459 


± 


0.140 


0.142 


19 


2 


Y 


-18.461 


± 


0.005 


0.273 


± 


0.160 


0.122 


13 


3 


J 


-18.552 


± 


0.002 


0.392 


± 


0.141 


0.123 


13 


3 


H 


-18.390 


± 


0.003 


0.302 


± 


0.109 


0.086 


13 


3 


Y 


-18.598 


± 


0.084 


0.844 


± 


0.418 


0.396 


11 


4 


J 


-18.665 


± 


0.092 


0.878 


± 


0.417 


0.384 


11 


4 


H 


-18.507 


± 


0.056 


0.808 


± 


0.276 


0.250 


11 


4 


Y 


-18.505 


± 


0.007 


0.139 


± 


0.471 


0.136 


7 


5 


J 


-18.606 


± 


0.003 


0.307 


± 


0.261 


0.125 


7 


5 


H 


-18.446 


± 


0.009 


0.538 


± 


0.273 


0.053 


7 


5 



Note. - Col. (1) Filter; Col. (2) Absolute magnitude for Am 15 (B) = 1.1 and 
no reddening; Col. (3) Luminosity vs. decline-rate slope; Col. (4) rms of fit in 
magnitudes; Col. (5) Number of SNe used in fits; Col. (6) Sample of SNe used in 
fit (see text). 
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Table 6. Fit Parameters 



Filter 


M x (01 


by 


Rv 


RMS 


Confidence level 


Subsample 




(mag) 


(mag decline-rate -1 ) 




(mag) 


(%) 




(11 


(21 


(31 


(41 


(51 


(61 


(71 


Y 


-18.418 ± 0.019 


0.495 ± 0.172 


1.7 


0.262 


99.6 


1 


Y 


— 18.535 ± 0.032 


0.611 ± 0.263 


3.1 


0.382 


98.4 


1 


Y 


-18.530 ± 0.000 


0.000 


3.1 


0.424 




1 


J 


-18.490 ± 0.020 


0.606 ± 0.170 


1.7 


0.258 


99.9 


1 


J 


-18.588 ± 0.028 


0.691 ± 0.230 


3.1 


0.336 


99.7 


1 


J 


-18.560 ± 0.000 


0.000 


3.1 


0.399 




1 


H 


-18.366 ± 0.014 


0.545 ± 0.131 


1.7 


0.201 


99.9 


1 


TT 

H 


— 18.432 ± 0.017 


0.581 ± 0.153 


3.1 


0.235 


99.9 


1 


H 


-18.389 ± 0.000 


0.000 


3.1 


0.292 




1 


Y 


-18.441 ± 0.005 


0.413 ± 0.146 


1.7 


0.145 


99.4 


2 


Y 


— 18.494 ± 0.007 


0.294 ± 0.187 


3.1 


0.184 


93.6 


2 


Y 


-18.503 ± 0.000 


0.000 


3.1 


0.189 




2 


J 


-18.514 ± 0.002 


0.376 ± 0.093 


1.7 


0.116 


99.9 


2 


J 


-18.574 ± 0.003 


0.363 ± 0.119 


3.1 


0.141 


99.6 


2 


J 


-18.545 ± 0.000 


0.000 


3.1 


0.163 




2 


H 


-18.390 ± 0.006 


0.469 ± 0.130 


1.7 


0.128 


99.9 


2 


H 


-18.415 ± 0.006 


0.459 ± 0.140 


3.1 


0.142 


99.8 


2 


H 


-18.398 ± 0.000 


0.000 


3.1 


0.169 




2 


Y 


-18.419 ± 0.005 


0.306 ± 0.115 


1.7 


0.098 


99.0 


3 


Y 


-18.461 ± 0.005 


0.273 ± 0.160 


3.1 


0.122 


94.1 


3 


Y 


-18.443 ± 0.000 


0.000 


3.1 


0.147 




3 


J 


-18.525 ± 0.003 


0.394 ± 0.113 


1.7 


0.101 


99.7 


3 


J 


-18.552 ± 0.002 


0.392 ± 0.141 


3.1 


0.123 


99.0 


3 
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Table 6 — Continued 



Filter 


M x (0) 


bx 


R v 


RMS 


Confidence level 


Subsample 




(mag) 


(mag decline-rate J 




(mag) 


(%) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


7 

J 




n nnn 

u.uuu 


3 1 


U. ±U^± 




Q 
O 


H 


-18.371 ± 0.004 


0.309 ± 0.097 


1.7 


0.073 


99.6 


3 


H 


-18.390 ± 0.003 


0.302 ± 0.109 


3.1 


0.086 


99.1 


3 


11 




n nnn 

u.uuu 


O. J- 


n 1 1 q 

U. 1 -Liy 




Q 
o 


Y 


-18.404 ± 0.063 


0.605 ± 0.323 


1.7 


0.302 


95.5 


4 


Y 


-18.598 ± 0.084 


0.844 ± 0.418 


3.1 


0.396 


96.2 


4 


V 

i 


1 o corv i n nrvn 

— lo.tjou zn u.uuu 


n nnn 

u.uuu 


^ 1 


n ^ns 

u.ouo 






J 


-18.494 ± 0.075 


0.746 ± 0.353 


1.7 


0.222 


96.7 


4 


J 


-18.665 ± 0.092 


0.878 ± 0.417 


3.1 


0.384 


96.7 


4 


T 

J 


i o C71 _(_ n nnn 

— lo.Oil It u.uuu 


n nnn 
u.uuu 


0. 1 


n m r 

U.OID 




^± 


H 


-18.339 ± 0.050 


0.591 ± 0.254 


1.7 


0.236 


97.6 


4 


H 


-18.507 ± 0.056 


0.808 ± 0.276 


3.1 


0.250 


99.1 


4 


u 
1 1 


i o q«q 4. n nnn 

— io.ouo in u.uuu 


n nnn 

u.uuu 


3 1 


n 






Y 


-18.464 ± 0.006 


0.004 ± 0.425 


1.7 


0.120 




5 


Y 


-18.505 ± 0.007 


0.139 ± 0.471 


3.1 


0.136 




5 


Y 


-18.507 ± 0.000 


0.000 


3.1 


0.155 




5 


J 


-18.591 ± 0.003 


0.245 ± 0.230 


1.7 


0.101 


83.9 


5 


J 


-18.606 ± 0.005 


0.307 ± 0.261 


3.1 


0.125 


85.8 


5 


J 


-18.628 ± 0.000 


0.000 


3.1 


0.143 




5 


H 


-18.428 ± 0.008 


0.444 ± 0.270 


1.7 


0.054 


91.5 


5 


H 


-18.446 ± 0.009 


0.538 ± 0.273 


3.1 


0.053 


94.9 


5 


H 


-18.442 ± 0.000 


0.000 


3.1 


0.106 




5 
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Note. - Col. (1) Filter; Col. (2) Absolute magnitude for Amis (B) = 1-1 an d no reddening; 
Col. (3) Luminosity vs. decline-rate slope; Col. (4) Total- to-selective absorption coefficient; Col. 
(5) rms of fit in magnitudes; Col. (6) Percent confidence level that slope is not zero; Col. (7) 
Sample. 
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Table 7. NIR Time of Maximum and Absolute Magnitude at Maximum 



r liter 


l B 


(M) 


a 


AT 

N S Ne 




(days) 


fmag) 

V o/ 


(mag) 

V o/ 




(1) 


(2) 


(3) 


(4) 


(5) 






Slow to mid-decliners 






Y 


-3.649 


-18.564 


0.129 


21 


J 


-3.178 


-18.651 


0.135 


21 


H 


-3.735 


-18.425 


0.130 


21 






Early-Peaking Fast-Decliners 






Y 


-2.500 


-18.471 


0.190 


1 


J 


-1.820 


-18.561 


0.200 


1 


H 


-2.140 


-18.328 


0.190 


1 






Late-Peaking Fast-Decliners 






Y 


4.230 


-17.668 


0.127 


2 


J 


2.490 


-17.669 


0.125 


2 


H 


3.950 


-17.797 


0.125 


2 



Note. - Col. (1) Filter; Col. (2) Time of NIR maximum; 
Col. (3) Average absolute peak magnitude; Col. (4) rms of fit in 
magnitudes; Col. (5) Number of SNe used in fits. 
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Table 8. Fit Parameters of Peak Magnitude Versus Ami^(B) for a Bimodal Distribution 



Filter X 


M x (0) 


bx 


RMS 


N S Ne 


Subsample 




(mag) 


(mag decline-rate -1 ) 


(mag) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Y 


-18.497 ± 0.010 


0.270 ± 0.160 


0.181 


20 


2 


J 


-18.551 ± 0.005 


0.332 ± 0.112 


0.146 


20 


2 


H 


-18.418 ± 0.008 


0.417 ± 0.123 


0.145 


20 


2 


Y 


-18.507 ± 0.022 


0.098 ± 0.261 


0.137 


8 


5 


J 


-18.617 ± 0.013 


0.240 ± 0.178 


0.124 


8 


5 


H 


-18.451 ± 0.012 


0.330 ± 0.173 


0.074 


8 


5 



Note. - Col. (1) Filter; Col. (2) Absolute magnitude for Ami 5 (B) = 1.1 and 
no reddening; Col. (3) Luminosity vs. decline-rate slope; Col. (4) rms of fit in 
magnitudes; Col. (5) Number of SNe used in fits; Col. (6) Sample of SNe used in 
fit (see text). 



-38 - 




20 40 
Time 



Fig. 1 — Left: Absolute magnitude B and V light-curves of SNe 2006D, 2004eo, and 2005el. 
The light-curves are shifted such that maxima agree. All SNe have similar decline-rates in 
the range of Amis(i?)=1.35 — 1.39 mag. Blue squares represent SN 2006D, green circles 
represent SN 2004eo, and red triangles represent SN 2005el. Note the similarity in the B and 
V light-curves. Right: Absolute magnitude YJH light-curves of the same SNe. Again, the 
light-curves are shifted to a common maximum. Note the difference in the YJH light-curves, 
especially the difference in strength of the second maxima in the J-band. 
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Fig. 2. — Left: Absolute magnitude B and V light-curves of SNe 2006mr, 2005ke, and 2007on 
(Am 15 (S) = 1.75-1.89). Blue squares represent SN 2006mr, green circles represent SN 2005ke, 
and red triangles represent SN 2007on. The light-curves are shifted to a common maximum. 
As in Fig. (U the optical light-curves are very similar. Right: Absolute YJH magnitude 
light-curves of the same SNe. The light-curves are shifted to a common magnitude. Note 
the large difference in the NIR light-curves, especially between SN 2007on, which shows a 
rise to second maxima. 
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Fig. 3. — Plot of the difference in peak magnitude derived from template vs. spline fits in 
YJH as a function of decline-rate for objects with observed NIR maxima. Uncertainties are 
smaller than the points, unless shown. Evidence of systematic differences are present, but 
the systematic errors are not large, compared with the uncertainties in the final absolute 
magnitudes. 
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Fig. 4. — Plot of the difference in peak magnitude derived from template vs. spline fits 
in YJH as a function of decline-rate for SNe whose NIR observations begin within 5 days 
after NIR maximum. Uncertainties are smaller than the points, unless shown. The weighed 
averages presented here are about the same as the weighted averages found in Fig. [31 
suggesting that when the observations begin within 5 days of NIR maximum, SNooPy does 
an adequate job of deriving peak light-curve parameters from template fits. 
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Fig. 5. — Fits of absolute YJH magnitude vs. decline-rate for subsample 1. Filled circles 
represent the SNe fit with a template, filled squares represent SNe fit with a spline. Open 
symbols represent fast-declining SNe and the highly reddened SNe. Uncertainties associ- 
ated with absolute magnitude are smaller than the points, unless shown. All three of the 
NIR bands show a statistically significant luminosity vs. decline-rate correlation, with a 
dependence in all three NIR bands comparable with those found in optical bands. 
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Fig. 6. — Fits of absolute YJH magnitude vs. decline-rate for subsample 2. The symbols for 
points are the same as those used in Fig. 5. A significantly weaker dependence of absolute 
peak magnitude on decline-rate is found for the YJH bands for subsample 2, and the rms 
dispersion decreasing significantly in all bands as well. 
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Fig. 7. — Fits of absolute YJH magnitude vs. decline-rate for subsample 3. The symbols for 
points are the same as those used in Fig. 5. Uncertainties associated with absolute magnitude 
are smaller than the points, unless shown. The l"-band shows a weak luminosity vs. decline- 
rate relation, while the J- and //-bands have a stronger relation, but the dependence of peak 
magnitude on decline-rate for subsample 3 is significantly less than in optical bands. 
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Fig. 8. — Fits of absolute YJH magnitude vs. decline-rate for subsample 4. The symbols 
for points are the same as those used in Fig. 5. This subsample produces the largest peak 
luminosity vs. decline-rate dependence of any set examined. 
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Fig. 9. — Fits of absolute YJH magnitude vs. decline-rate for subsample 5. The symbols 
for points are the same as those used in Fig. 5. Uncertainties associated with absolute 
magnitude are smaller than the points, unless shown. There appears to be little dependence 
of peak luminosity on decline-rate within the uncertainties for the Y- and J-bands and 
marginal dependence for the .ff-band. These results are consistent with the luminosity vs. 
decline-rate relations derived from subsamples 2 and 3. 
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Fig. 10. — Top: Hubble diagram constructed from SNe and fit parameters of subsample 
2. The blue circles represent the averaged NIR distance moduli for each SN, and the red 
triangles represent the i?-band distace moduli for each SN. The solid line shows the standard 
cosmology redshift-distance relationship (eq. [2]). Middle: The NIR residuals with respect to 
the standard cosmology model. Bottom: The 5-band residuals with respect to the standard 
cosmology model. The decrease in scatter from the optical bands to NIR bands provides 
strong evidence that SNe la in NIR bands are excellent standard candles. 



